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The stabihty of the flux hne lattice has been investigated within anisotropic 

(N . 

^ . London theory. This is the first fuh-scale investigation of instabihties in the 

. 'chain' state, the equilibrium lattice that is similar to the Abrikosov lattice 

«) , at large fields but crosses over smoothly to a pinstripe structure at low fields. 

0> , By calculating the normal modes of the elasticity matrix, it has been found 



the lattice is stable at large fields, but that instabilities occur as the field is 
reduced. The field at which these instabilities first arise, b*{€,9), depends 



c . 

O . on the anisotropy e and the angle 9 at which the lattice is tilted away from 

the c-axis. These instabilities initially occur at wavevector k*(e, 0). The 
dependence of k* on e and 6 is complicated, but the component of k* along the 
average direction of the fiux lines, kz, is always finite. For rigid straight flux 
lines, the cutoff necessary for London theory has been 'derived' from Landau- 
Ginzburg theory, where the shape of the vortex core is known. However, for 
investigating instability at flnite kz it is necessary to know the dependence of 
the cutoff on kz, and we have used a cutoff suggested by Sudb0 and Brandt. 
The instabilities only occur for values of the anisotropy e appropriate to a 
material like BSCCO, and not for anisotropics more appropriate to YBCO. 
The lower critical held Hc^{(j)) is calculated as a function of the angle (p at 
which the applied field is tilted away from the crystal axis. The presence of 



kinks in Hc^{4>) is seen to be related to instabilities in the equilibrium flux 

line structure. 

PACS: 74.60.Ec, 74.60.Ge 
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I. INTRODUCTION 



The existence of novel flux line structures in high-temperature superconductors has led 
to intensive investigation of the mixed-state of these materials. Unusual structures have 
been observed in Bitter pattern experiments on YBCO0I and BSCCOii when the apphed 
magnetic field was tilted away from the c-axis. In YBCO, the 'chain' state was observedBi 
where the anisotropy of the material causes the usual repulsive flux line interaction to become 
attractive within the tilt plane, the plane containing the magnetic field and the c-axis. The 
chain state had been predicted within the framework of the London approximationii. In 
some of the experiments! the chains of flux lines were seen embedded in an approximately 
triangular flux line lattice, but it is believed the presence of the lattice was due to pinning of 
the flux lines. In BSCCO, similar structures of chains embedded in a lattice were also seeniS, 
but the dependence of the flux line spacings on the tilt angle and magnetic field were different 
from those in YBCO, implying these structures may be created by a different mechanism. 
Possible explanations for the flux line structure seen in BSCCO have been proposedSi. It 
was suggested there existed inter-penetrating flux line lattices, one orientated approximately 
parallel to the c-axis while the other is orientated approximately parallel to the ab-plane. 
Within the framework of the London approximation, it has been shown that provided the 
anisotropy is large enough there is more than one possible angle at which the flux lines 
initially enter the samplei. 

The mixed state for isotropic superconductors is a periodic triangular array of straight 
flux line^. In uniaxially anisotropic superconductors it was predicted0'0 the flux lines 
would form a distorted triangular lattice, where the spacings between the flux lines depend 
on the strength of the magnetic field B and the anisotropy mass ratio M^jM. A lattice 
similar to this distorted triangular lattice has been observed in YBCO at large fields using 
small angle neutron scatteringS. 

The stability of this distorted lattice against elastic deformations has been studied exten- 
sively within London theory. Sudb0 and Brandt@ showed that at large anisotropy M^/M 3> 
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1 and small magnetic induction b = B/Hc2 ^ 1 theenergy associated with a pure shear- 
ing mode of the flux lattice canbecome negative. The existence of a tilt-wave instability, 
k = (0,0,^2), was demonstrated by Sardella and Moore0 and confirmed by Nguyen and 
Sudb0 0, who both employed the =ame cutoff procedure. The distorted trtogular lattice is 
the lattice one would expect using the ideas of anisotropic scalingEl. The lattice's basis vec- 
tors are proportional to 1/ \/h and depend in a simple manner on 9, k, e. Minimizing the free 
energy, Daemen et al.0 showed that out of the set of centered rectangle lattice structures, 
one of which is the distorted triangular lattice, the true equilibrium lattice behaves quite 
differently. At large fields it is approximately the distorted triangular (Abrikosov) lattice, 
but there is a smooth crossover to the 'chain' state at small fields. This state has one of the 
basis vectors independent of the field, with the other being inversely proportional to h. 

In this paper we investigate the stability of this equilibrium lattice at a general wavevector 
k = {kx, ky, kz). The existence of a zone center instability has been observedlll, but this is 
the first full-scale investigation of elastic instabilities of this equilibrium lattice. The cutoff 
used is that proposed by Sudb0 and BrandtEl, which depends on kz- The lattice is found 
to be stable at large fields. As the field is reduced, the field at which the instability first 
occurs b*{e,6) depends on the anisotropy e, and the angle 6 at which the lattice is tilted 
from the c-axis. The magnitude of these fields b = B / H^^ ~ O(10~^) is approximately that 
used in the experiments where the unusual flux line structures were seen. It is generally 
believed these instabilities are somehow related to the different flux line structures seen in 
experiments, but we have as yet no good theoretical interpretation of the connection. 

The attractive vortex interaction also forces the flux lines to initially enter a supercon- 
ductor not as single flux lines but in chains. For very large anisotropy, there exists the 
possibility that these chains may first enter the sample at more than one direction. The 
precise details are complicated, see section [V^. However, they indicate that for YBCO the 
chain state is the stable low field structure. This is not the case for BSCCO, where it is 
possible the low field structure will be the coexistence of different chain state orientations. 

The initial instabilities observed always have finite k^. For instabilities of non-zero kz 



the form of the cutoff used within the London theory is cruciaL In the calculation of Nguyen 
and Sudb0 111, which used a cutoff that did not depend on kz, the coexistence of different 
flux line orientations and the zero field tilt wave instability were two different effects, i.e. 
they initially occurred at different anisotropies. In the limit 6 ^ it can be shown that 
the onset of both effects are related to the line tension Pi{6) = si{9) + d'^Si/dO'^ becoming 
negative. If the cutoff that depends on kz is used, it is found that the two effects do occur 
at the same anisotropy, implying this cutoff may be more reliable. 

II. LONDON THEORY 

A convenient way to describe the low- field magnetic properties of high- Tc superconductors 
is London theory. In the isotropic form, this theory just depends on the penetration depth 
A and on k = A/^, where ^ is the coherence length. To allow for the anisotropy of the HTSC 
compounds, the square penetration depth A^ is replaced by the tensor A^^ = Ajj. Here, we 
shall only investigate uniaxial anisotropy where Kxx = Ayy = A^^ 7^ Azz = A^ are the only 
non-zero elements of Aij. The anisotropy of the material is governed by the parameter e 
which in the effective mass model is given by = Mxx/Mz, and \ab/K = idiab = e- 

The London free energy can be written as 




where A is the vector potential of the magnetic field H = V x A, (^9 is the phase of the 
order parameter, $0 is the flux quantum, and A~^ is the inverse of the square penetration 
depth tensor A. The magnetic induction B is the average magnetic field B = (H) = i?z. 
London theory is a good approximation at low induction, B < Q.2Hc^, where the cores do 
not overlap strongly. 

In general, the B field is not aligned with the crystal axis, and we chose our coordinate 
system such that the B field lies in the X — Z plane and is tilted away from the c-axis by an 
angle 9. It is then convenient to use the 'vortex' coordinate system (xyz). This is obtained 
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by rotating the crystal frame (XYZ) by an angle 9 around the Y-axis, see Fig. 1. In the 
vortex coordinate system, the square penetration depth is given by A^^ = Ai^^/j + K-iCaCf^ 
where Ai = A^^, A2 = — A^j,, (a,/?) = (x,?/, z), and Cq is the a component of the unit 
vector c in the vortex frame. 

The free energy can be written in a simpler form. Minimizing ([^) with respect to the 
vector potential A, and then taking the curl of the equation, we obtain the London equation 

H + Vx{A.VxH} = $o5: 1 dr,53 (r - r,) (2) 

i 

It is possible to derive this equation from the Ginzburg-Landau equations, assuming the 
order parameter has constant magnitude. The right hand side of the London equation 
comes from V x Vy^ = I]j / drj^a (r — rj), where rj is the position of the i flux line. 
The London equation is linear, and has the solution 

H„ (r) = %Y.J (r - r,) (3) 

i 

Using this potential, the London free energy ([^) may now be written as 

F=^T.f / dr,"drf (r, - r,) (4) 



The convention of summation over repeated indices is assumed. Written in this form, the 
free energy can be seen as consisting of two parts, the self-energy terms { i = j ) and the 
interaction terms {i 7^ j). 

The Fourier transform of the potential V{r — r^) is0 

A2gag/3 



Vaf, (k) 



1 + AiP 



Sa(3 — 



(5) 



1 + AiA;2 + A2g2 

where q = k x c. From Eq. (^), we see that the potential decays only as l/A;^ as A; 00, 
implying that H (r) is singular at r = r^. The singularites in H (r) are due to the absence 
of the vortex cores from London theory. A convenient way to circumvent the problems 
associated with the divergences is to introduce a cutoff into the London potential (|^) via 



T/ (^.^ ^^^^ 

Vaf3 (k) = - 



AiF 



1 + AiA;2 + A2g2 



(6) 



This is equivalent to replacing the delta function in the London equation with a short- 



ranged function S {r — rj). By making this short-ranged function a gaussian, of width y/2^ab 
along a and b and width ^/2^ along c, the cutoff becomesEl'S 

S{k) =exp-2g(k) 
g{k) = C\ (k X c)^ + C (k ■ c)^ (7) 

= C^^g' + Cik' - q') 

This provides an elliptical cutoff at large k^ = {kx, ky), as required from the shape of the core 
in Ginzburg-Landau theoryS. The factor 2 in the exponential of 5'(k) is just convention, 
but it can be determined more accurately by comparison with results from Ginzburg-Landau 
theory. The results in this paper are not affected by the choice of this parameter. 

The exact form of the cutoff is of some debate. For straight rigid vortices the form 
of the cutoff can be derived from the shape of the core within Ginzburg-Landau theory^. 
However, the cutoff in Eq. (|^) depends on k^, the component of k in direction of the B 
field. Sardella and Moored and Nguyen and Sudb00 investigated the tilt-wave instabilities 
of the distorted triangular lattice using a cutoff that depended only on k_|_ = {kx,ky), i.e. 
g(k) = g(k^) = ^l^ik^ x c)^ + ^^(k^ ■ c)^ in Eq. (g). We believe the cutoff in Eq. (g) is 
more physical, as it does not depend on being able to specify k^. That can only be done by 
reference to the average direction of the flux lines, i.e. B, but it is hard to believe that the 
cutoff should be sensitive to this overall average direction. Some authors do not use a cutoff 
function ^(k), but instead introduce an upper limit on any integrations over k. In most 
situations, if the symmetry of the upper limit introduced is the same as the cutoff function 
S(k), then similar results are obtained^. However, in some situations care may be required 
to ensure it is does not matter whether certain points in k-space are just inside or outside 
the integration ranged. 
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III. EQUILIBRIUM LATTICE 



In an isotropic superconductor, the equilibrium flux line lattice is a periodic array where 
the unit cell is defined by an equilateral triangle. Defining the x-axis to coincide with one 
of the basis vectors, the basis vectors may be written as Ri = ax and R2 = a{± + V3y)/2 
where = 2<^o/VSB. The presence of anisotropy causes the lattice to distort from this 
equilateral structure. Using the ideas of anisotropic rescalin^li, the lattice is expected to be 
of the form 

Ri = a7x 

R2 = a(7x + y3y/7)/2 (8) 

where 7"^ = cos^ 9 + sin^ 9. This structure is indeed the equilibrium structure in the limit 
of the lowest Landau level, and was seen to minimize the free energy of a set of rescaled 
structuresii. All the length scales depend on the strength of the magnetic field in a similar 
manner i.e. they are proportional to l/Vb. We refer to this rescaled structure as the 
Abrikosov lattice. 

The presence of anisotropy dramatically changes the profile of the magnetic field. If the 
anisotropy is large, the magnetic field associated with a single isolated fiux line contains 
regions around the fiux line where the local magnetic field points in the opposite direction 
to the average field. This allows the usually repulsive fiux line interaction to be attractive. 

Daemen et al.0 showed if we investigate the set of fiux line structures with a centered 
rectangular symmetry, one of which will be the same structure as the Abrikosov lattice, 
different equilibrium structures can be seen. The unit cell consists of an isosceles triangle, 
see Fig. 2, with two sides of length I2 the other length li, with an angle ip between sides of 
length li and I2. The lattice vectors are 

Rmn = i^h + COS 1p)^ + ^^2 siu tpf (9) 

where m and n are integers. Repeated computer minimizations of the total free energy have 
confirmed that for uniaxial anisotropy the assumption the unit cell is an isosceles triangle is 
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validS. The magnetic flux per unit cell must be one flux quantum, which allows /i, I2 and 
ijj to be written in terms of one parameter 



h 
h 

cos?/' 



$0 



p 



h/p 

1 



1/2 



(10) 



Following Daemen et al.0 we minimize the free energy per unit cell e using the golden- 

because the derivatives of e with respect to p are hard to calculate. 

The dependence of li and I2 on the field, h = BjHc^iO = 0), is shown in Fig. 3. At 
large fields the two lengths scale as approximately 1 / ^/b but there is a smooth crossover to a 
regime at lower fields where one is approximately constant while the other scales as 1/6. The 
region where this crossover occurs is characterized by n, 6 and e. The low field equilibrium 
state corresponds to the chain state, or 'pinstripe structure' observed by Gammel et al.0 in 
YBCO. This shows care must be taken in defining the flux line lattice in any calculation, 
as in this low field regime the Abrikosov state is not an equilibrium state, see Fig. 4. The 
value of p corresponding to the anisotropically rescaled Abrikosov lattice is given by pabr- 

The reciprocal lattice of this equilibrium lattice has basis vectors 

1 . 



_ 27r. 



■ 271 2n 

m- i^~r~ cosip 

h h 



sin-?/'' 



(11) 



In the calculations that follow, wavevectors will be measured in units of 
(27r/ (3/i)x, tt/ {I2 sin ip)y), making the rescaled reciprocal lattice vectors Qmn = 3rix+ (2m — 
n)y. 



IV. ELASTIC THEORY OF FLUX LINE LATTICE 

The minimum free energy configuration of the flux lines has been assumed to be a periodic 
array, whose unit cell is an isosceles triangle with the base orientated along the x-axis. To 
check the assumption that this free energy is at least a local minimum, the change in the 
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free energy associated with small displacements Sa (Ri{z)), a = {x,y), from the equilibrium 
lattice Rj = nUi + mR2 can be derivediH. Keeping terms only to second order in the 
displacements SQ,(Ri) the change in the free energy is 

AF = ^ 1 d'ks^ (-k) ^^p{k)sp{k) (12) 

where the integration over k_L = {kx, ky) runs over the first Brillouin zone and over kz on 
the interval (— cxd,cx3). The elasticity matrix is 

(k) = E iU (k + Q) - U m (13) 

Q 

where 

/a/3 (P) = P^K/3(P) +PaPpVzz{v) - PzPcMzM 

-PzPpVzM (14) 

The stability of a periodic lattice may be determined by investigating whether the normal 
modes of the elasticity matrix ^^^^(k) always remain stable i.e. the eigenvalues of $Q,^(k) 
are positive, or whether in some regions of the Brillouin zone the normal modes become 
unstable. 

The stability of the distorted triangular lattice (Abrikosov lattice) has been examined 
by various authorsEHll. Sudb0 and BrandtEl observed that for a configuration of rigid flux 
lines, i.e. kz = 0, as the magnetic field was reduced below a specific level the normal modes 
became unstable. Sardella and Moore0 observed a tilt wave instability, k^ = ky = 0, kz 7^ 0, 
and this instability was present in all fields. However, the cutoff used did not depend on kz- 
It was suggested that with the use of the cutoff of Eq. (0) the instability would disappearii. 
For k = (0, 0, kz) the elasticity matrix is diagonal, and the instability calculated by Sardella 
and Moore was associated with becoming negative. Upon repeating the calculation of 
Sardella and Moore with the cutoff (|^), we found that the eigenvalue ^yy did indeed remain 
positive. However, became more unstable, and as in Ref.lii, the lattice becomes unstable 
as the field is reduced beyond a critical level. Nguyen and Sudbji have confirmed this by 
calculating the normal modes of rigid flux lines in the limit 6 = 0. 
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The Abrikosov lattice is only a good approximation of the equilibrium structure at fields 
b > O(10~^ — 10^^). We have investigated the normal modes of the elasticity matrix $o^(k) 
for k = {kx, ky, kz), with the equilibrium lattice being the isosceles triangle described in 
Section |TI|, using the kz dependent cutoff (^. The presence of instabilities is expected as 
SardellaS has confirmed the existence of a zone center instability for large tilt angles 9, large 
anisotropy 1/e and small fields b. 

The fields at which instabilities are observed depend on k, on the tilt angle 6 and on 
the anisotropy e. The isotropic system, e = 1, is stable. As the anisotropy increases, e 
decreases, no instabilities are observed until a critical anisotropy is obtained. The minimum 
anisotropy required for instabilities to be seen depends on k. For k, = 20 instabilities were 
seen for > 120.3, but this was increased to 1/e^ > 138 for k = 50. For values of 1/e less 
than these critical values the lattice is stable at all orientations. 

If the anisotropy is larger that this critical anisotropy, instabilities can be seen. The 
lattice is always stable at large fields. As the field is reduced there is a specific field b*{n, e, 9) 
at which the lattice initially becomes unstable. Fig. 5 shows that for a given angle 6 there 
is a minimum anisotropy 1/e below which the lattice is stable in all fields. It also shows 
that the critical field b* increases as the system becomes more anisotropic. The lattice is 
always stable at ^ = and 6 = 7r/2, but once the anisotropy is large enough there exists a 
range of angles at which instabilities are seen, see Fig. 6. This range of angle increases as 
the anisotropy increases. 

The wavevectors at which these instabilities are first seen, k*(/t, e, ^) always have finite 
kz- Both b* and k* are functions of k, e and 6. These unstable modes correspond to 
displacements of the lattice approximately parallel to the x-axis, and we call it the 'staircase 
wave' instability to distinguish it from other instabilities observed. In Fig. 7 the dependence 
of the X and y components of k* are shown; k^ and ky have been rescaled so that the dashed 
line shows the edge of the first Brillouin zone. 6 decreases from small k^ to large kx- The 
explanation of the behavior of the actual wavevector where the instability first appears, k*, 
is not obvious to us. 
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The parameters chosen in this paper (e = 1/60, k = 50) were chosen to give some insight 
into the behavior of BSCCO. While the exact value of the anisotropy in BSCCO is unknown, 



the values used in this paper are similar to those used in other paper^i20J13. In comparison. 
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YBCO is much less anisotropic, with e ~ 1/5 and k ^ 50, and as can be seen from Fig. 5 is 
well below that required for the instabilities to be observed. To explore the instabilities in 
full detail as a function of k, e, 6 and b is very time consuming, but Fig. 5 shows the 'chain' 
state is stable for parameters that could describe YBCO, but unstable for BSCCO. 

V. LOWER CRITICAL FIELD: ABRIKOSOV LATTICE LIMIT 

The presence of elastic instabilities for very anisotropic materials may indicate why there 
are differences in the Bitter patterns observed for BSCCO and YBCO. One interpretation 
of the BSCCO patterns was the coexistence of two interpenetrating latticesi'l. This pos- 
tulate has been investigated by calculating the angles at which single flux lines first enter 
a superconducting sample^lS. We modify the calculation of Nguyen and Sudb(70 to show 
this possible effect in order to emphasize the importance of the cutoff procedure used. 

We consider a cylindrical superconducting sample, with the applied field H perpendicular 
to the axis of the cylinder, and tilted at an angle away from the crystal c-axis, see Fig. 
8. This geometry is chosen so that demagnetization effects permit solutions where the flux 
lines are straight for all orientations of the applied field. 

Within this geometry, the Gibbs free energy for a system of rigid straight flux lines, tilted 
at (a different) angle 9 is 

^^P_B2UMl^ (15) 

Neglecting the interaction between the flux lines, the free energy within the London approx- 
imation (^) is 
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where Xg = A^^ sin^ ^ + -^c cos^ 9. As with most calculations within the London approxima- 
tion, the integral is formally divergent without the cutoff term, but the use of the cutoff ^(q) 
described in the previous section allows the calculation to proceed. Also, as the integration 
is over q perpendicular to the flux line this calculation is insensitive to the dependence of 
the cutoff S'(q). By neglecting the interaction between the flux lines we are essentially just 
discussing the behavior of a single flux line. This calculation can be viewed as the extension 
of the Abrikosov lattice to 6 = 0. 

As the field Ti is increased, the flux lines initially enter the sample when G = 0. Although 
it would be preferable to chose 0, the orientation of the applied field, and then calculate 9 
and H it is only possible to assume 9 and then calculate the corresponding values of H and 

0. Looking for solutions where dG/d9 = implies the orientation(s) of the flux lines within 
the sample are governed by the relation 

tan0 = ^ (17) 

^ 1-F'tan9/F ^ ' 

When investigating the first entry of flux lines into the sample, it is well known that the 
cutoff used is important. If a circular cutoff is used then there is a unique orientation of the 
flux lines, given byiS 

tan0 = e^tan6' (18) 

However, this is not true if an elliptical cutoff is used, as for large anisotropies there exists 
a nonmonotonic relation between 9 and 0, see Fig. 9. This implies the possibility of the 
existence of two orientations of flux lines. 

Within the nonmonotonic regime, there exist three possible orientations of the flux lines, 
6'a, 9h, and 9c, for any given 0. 9b corresponds to an unstable orientation of the flux lines 

1. e. a maximum in the Gibbs free energy, and will not be considered further. 

By calculating the Gibbs free energy, we can see whether the flux lines will orientate 
at either 9a or 9^, and the relationship between and 9 determined. A 'forbidden' region 
occurs due to the presence of the nonmonotonic regime. The flux line lattice can only be 
orientated at angles 9 < 91 and 9 > 92. 
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This can be seen by noting the relationship between the calculated values of Ti. and 
(f). Fig. 10 shows 7i(0) when all flux line orientations are allowed. The points in Fig. 
10 are equally spaced in 6. For large anisotropy 7i(0) becomes nonmonotonic when 6{(j)) 
becomes nonmonotonic. In Fig. 11 the three values of 7i(0) in the nonmonotonic regime 
correspond to the three possible orientations of the flux lines at 6a, 6b and 6c. While all 
three orientations correspond to G = and dG/d6 = 0, only one is a global minimum in 
G{6). This global minimum will be the equilibrium structure and is the orientation with the 
smallest value of The lower critical field Hc^{(f)) is therefore just the smallest value of 

7i(0) for any given (f). From Fig. 11 it can be seen Hc^ has a kink where 1-L{6\) = 7i(^2) ^^^i 
4>{6l) = 0(^2) ~ (pkink- The presence of kinks in the lower critical field for single flux lines 
have also been observed within the Ginzburg-Landau model0. 

If we now allow the superposition of non-interacting flux lines, which are assumed to 
be far apart in the limit 6^0, the 'forbidden region' may be removed. At Hc^{(j)kink) the 
average field B may be orientated at all angles within the 'forbidden' region 6^ < 6 < 62 hj 
orientating some flux lines at 6^ and the others at 62- 

Various authors have discussed the competition between the coexistence of different flux 
line species and the presence of elastic instabihties in the supposed equilibrium structure. 
For single isolated flux lines (in the limit 6-^0) the elastic instabilities depend most strongly 
on kz and only tilt-wave instabilities will be discussed further in this section. 

As kz is non-zero, the cutoff used is important. Nguyen and Sudb0 used S'(k) = ^(kx). 
For any given value of k, they found two important values of the anisotropy. At small 
anisotropics the relationship between and 6 is monotonic. For 1/e > T^^ the relationship 
between and 6 becomes nonmonotonic and the Gibbs free energy may be doubly degen- 
erate. However, if 1/e > F^"^ the 'perpendicular' tilt modulus c^4^(k) = ^xx{0,0,kz)/kl 
becomes negative over a range of angles. The two critical anisotropics depend linearly on 
luK, but r^'^ is much larger than T^^. 

There is a marked difference if the cutoff >S'(k) = exp{—2g{'k)), gCk) = ^^^(k x c)^ -|- 
(k ■ c)^. The dependence of T^^ and F^"^ on k, is similar, but F^"^ = F^*^ = F. For 
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K = 50 and 5'(k) = exp{—2g(k)) it is found F ^ 9.65. The value of T does depend on 
the cutoff procedure used. Using a cutoff of the same symmetry but a different strength, 
Sik) = exp(-^(k)), it is found T ^ 10.03 for k = 50 

These two 'competing' effects are both related to the line tension Pi{6) = ei{6) + 
d'^ei{6)/d6'^ , where e is the line energy. In the limit 6-^0 where the London free energy 
F = nei, n is the areal density of flux lines, then it follows from ( p!7|) that 

P,{9)=eii9)sec'i^-9)^ (19) 

However, C44 may be defined as C44 = (PF{9)/d9'^. For the isolated noninteracting flux 
lines F = nei, so C44 oc ei{9) + d'^ei{9)/d9'^ = Pi{9)^, where the first term is due to the 
compression of the flux line lattice during tilting. Therefore the presence of a kink in Hc^ (0) 
not only causes a restriction of orientations of the flux line lattice, and the possibility of the 
coexistence of different flux line orientations, but is also related to the presence of a tilt- wave 
instability. The tilt-wave instabilities only occur where d(j)/d9 < 0, see also Grishin et al.i 
, but this is only part of the region of excluded flux line orientations 91 < 9 < 92 and the 
instabilities are never observed. 



VI. LOWER CRITICAL FIELD: CHAIN STATE LIMIT 

The previous calculation assumed the equilibrium low field flux line structure was a 
configuration of well separated non-interacting flux lines ,i.e. the extension of the Abrikosov 
lattice to 6 = 0. In section |T| it was seen the low fleld equilibrium structures are very 
different from the Abrikosov lattice. The anisotropy induced attractive interaction makes it 
easier for chains of flux lines to enter a sample than for single flux lines, which forces ifci for 
an single inflnite chain to be less than ifci for a single flux linJ. We take the limit 6 — > to 
imply there is only one vortex chain in the sample, the separation between the chains being 
effectively inflnite. 

The presence of instabilities in the 'chain' state were investigated in Section which 
seemed to indicate a different behavior for materials such as YBCO from those such as 
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BSCCO. These instabilities were seen as the field was reduced from a stable regime similar 
to the Abrikosov lattice. In this section we investigate what happens as the fiux lines initially 
enter the sample i.e. b = 0. Chosing the same geometry as section |V] we again look at the 
relationship between (p and 6 and relate this to the elastic instabilities. 

For a single fiux line in section |V|, the 2D line energy integral is easily calculated by 
rescaling the coordinates, see Nguyen and Sudb0 0. For the single infinite chain, this 2D 
integral is replaced by a sum over the reciprocal lattice vectors of the Id chain and a Id 
integral perpendicular to the chain. Assuming the chain to be a set of fiux lines equally 
spaced along the a;-axis, then 

J (27r)2 U^J 2n ^ ^ 

where = 27m/ Ich, n = ± 1 ± 2... , Ich being the separation of the fiux lines in the chain. 
Ich is determined by finding the minimum in the line energy, and is a nonmonotonic function 
of ^. 

The behavior of the chain state is subtly different from that of the single flux line. The 
main effects are the same, i.e. for large anisotropy there exist a range of angles over which 
the chains cannot be orientated, but the details are more complicated. 

The chain state is stable for small anisotropics, but for large anisotropics the relationship 
between and 6 once again becomes nonmonotonic. For this to occur it requires a larger 
anisotropy than for the single flux line. Fig. 12 shows (j){6) just after the onset of nonmono- 
tonicity for isolated flux lines. The chain state has no restriction on the orientation of the 
flux lines, unlike the single flux line, and also the chain state always has a lower (or equal) 
lower critical field Hc^, see Fig. 13. 

Increasing the anisotropy. Fig. 14, the chain state eventually shows signs of instabilities. 
Initially (j){6) becomes nonmonotonic at large tilt angles, 6 ~ 97r/20, and the chains are 
excluded over the range < 9 < 02. This occurs for 1/e > Ti{k) and Hci{4>) contains a 
kink in a manner similar to the single fiux line. 

However, for 1/e > r2(/t), 0(^) is also nonmonotonic at smaller angles, 6 ~ 7r/4. This 
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implies the chains are also excluded over a different range 9^ < 9 < Oy. Prom TC{4>), Fig. 
15, it can be seen that a stable range 91, < 9 < O^, exists between these two forbidden 
regions. The lower critical field Hc^[(j)) now has two kinks. As the anisotropy is increased 
further these excluded regions grow until at 1/e > r3(«;) the chain state is excluded over the 
whole range 9\ < 9 < 9^. The second kink in Hci{4>) disappears, Fig. 16, and only one kink 
remains. For n = 50, Fi a; 15, F2 ~ 30 and F3 a; 60. 

The elastic instabihties again only occur when d(j)/d9 < 0. These instabilities have 
slightly different properties at small 9 and large 9. At large 9 the lattice is most unstable 
to a tilt-wave instability k = (0, 0, kz) while it is a staircase wave instability, with kx and 
kz both non-zero, that causes the chain to become unstable at smaller 9. Again, these 
instabilities only occur in the forbidden regions, and the chain state is stable at all allowed 
orientations. 

VII. CONCLUSION 

We have investigated the presence of instabilities in the flux line lattice in anisotropic 
superconductors using London theory. There are limitations on the applicability of the 
London theory, but the wavevectors at which the instabilities appear are well within its 
limits of validity, i.e. k <^ 0(1/^). 

Due to the anisotropy induced attractive flux line interaction, the form of the equilibrium 
flux line lattice has to be determined numerically at all flelds. There is a smooth crossover 
from a large fleld Abrikosov-like state to the low fleld chain-like state. The nature of the 
instabilities observed in this equilibrium lattice is not complex. At large fields it is stable. 
As the field is reduced, at a particular field b* some of the normal modes of the elasticity 
matrix become unstable. This field depends on the anisotropy e and the angle 9 at which 
the fiux lines are orientated to the crystal c-axis. The instability is characterized by always 
having finite k^ and k^, and we call this a staircase wave instability to distinguish it from a 
tilt wave instability which is defined as an instability which depends only on k^. There is a 
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minimum anisotropy required before these instabilities are observed. This depends on k and 
for K = 20 the instabihties were seen for > 120.3, but this was increased to 1/e^ > 138 
for K = 50. This may indicate why the Bitter patterns observed on YBCO are different from 
those on BSCCO. 

As the component of k along the flux lines is non-zero, the choice of cutoff is crucial. 
We have used the cutoff suggested by Sudb0 and Brandt, where 5* (k) = exp— 25f(k) and 
g{k) = (k X c)^ + (k ■ c)^. This cutoff appears more physical than a cutoff that just 
depends on k^ and allows the elastic instabilities and the nonmonotonic behavior of (j){6) 
to be related to a single quantity, the line tension Pi{0), as expected. However, it should 
be noted there are still problems associated with this cutoff in the limit oo. Once 

in the region where kz forces 5* (k) -C 1, the elasticity matrix (|13]) will be dominated by 
— J^qfaisiQ)- While is zero, the eigenvalues of the elasticity matrix are negative, 
showing London theory is always unstable at all fields and angles. However, these instabilities 
occur at values of kz ~ 0(1/^) and are unphysical. The inclusion of core effects, e.g. core 
bending energy, may remove remove this instability, but this is outside the domain of validity 
of the London approximation. 

Instabilities can also be observed by investigating the lower critical field ifci, i-G. b = 0. 
These 6 = instabilities are different from the large field staircase wave instabilities. When 
the applied field 7i is tilted away from the crystal axis, the B-field and the applied field Ti. 
are not parallel. In Section it was seen that for small fields and large anisotropies there 
existed a nonmonotonic relationship between the angle 6 at which the B-field is tilted away 
from the c-axis, and the angle at which the applied field is tilted. The nonmonotonicity 
is related to the elastic instabilities observed, but also imply there is a restriction on the 
allowed orientations of the B-field and a kink in Hc^{(f)). The properties of the chain state 
changes at three different values of the anisotropy. For 1/e < Fi, all possible orientations 
of the chain state are possible, while in the interval Fi < 1/e < F2 the chains cannot be 
orientated over the range < 6 < 62. However, if F2 < 1/e < F3 the flux lines are excluded 
over two separate regions which grow to one large region for 1/e > F3. This unusual behavior 
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can also be described as Hc^{(j)) developing one kink at 1/e = Fi, a second kink developing 
at 1/e = r2 but Hc^{(f)) only having a single kink for 1/e > T^. 

Whether the peculiar effect of having the allowed orientations 9y<9< between 
two forbidden regions really exists is unclear. Just as the existence of the forbidden region 
9\ < 9 < 62 for the isolated flux lines in Section |V| may have been an indicator of a state 
with a lower free energy and applied field ifci, there may exist a new composition of flux 
lines that will be stable at large anisotropics and have a lower free energy and applied field 
i^d than the chain state. This may be a completely new state or a superposition of chain 
states and other flux lines, but a full investigation of such states is left for future work. 
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FIG. 1. The 'vortex' frame {x,y,z) is obtained by rotating the crystal frame {X,Y,Z) by an 
angle 9 about the Y axis. 




22 









7t/4 — 




57T/8 






I1 



10 ^ 10 10 ^ 10 ^ 10 

b 



FIG. 3. The field dependence of the separation of the flux lines (measured in Xab) for e = 1/60 
and K = 50. At large fields h and I2 are proportional to l/\/6 but at lower fields there is a crossover 
to the 'chain' state. 
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FIG. 4. The dependence of the energy per flux line on the parameter p, for e = 1/60, ^ = Stt/S, 
K = 50, b = 10~^. The arrow marks the value of p of the rescaled Abrikosov lattice. This figure 
clearly shows that the Abrikosov lattice is not the solution of minimum free energy. 
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FIG. 5. The field at which the lattice initially becomes unstable as a function of the anisotropy 
e, for 9 = 37r/8 and k = 50. 
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FIG. 6. The field at which the lattice initially becomes unstable as a function of the angle, 0, 
by which the lattice is tilted from the c-axis, for different values of e and k, = 50. 
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FIG. 7. The rescaled kx and ky components of k*, for different e and k = 50. The dashed Une 
marks the edge of the first Brillouin zone. 
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FIG. 8. The geometry of sample and appUed field H. The applied field is tilted at an angle (j) 
from the c-axis, while the magnetic induction B is tilted at an angle 9. 



25 





FIG. 10. The applied field 7^(0) for e = 1/5 and n = 50. 
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FIG. 11. The applied field for e = 1/60 and k = 50. 




FIG. 12. The relationship between ^ and 9 for different lattices, with k = 50 and e = 1/10. 
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FIG. 15. The relationship between and ^ showing the existence of two kinks in Hci{4>), for 
K = 50 and e = 1/40. This allows the orientation of the chains state between two forbidden regions. 




0.005 0.01 0.015 0.02 

(p/n 



FIG. 16. 'H{(p) for At = 50 and 1/e = 100, showing that at large anisotropy the second kink in 
i^ci disappears. 
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